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Abstract 

We investigate the periodic Anderson-Holstein model by using the dynamical mean-field theory combined with the 
exact diagonalization method. For the strong electron-phonon coupling g Z gc, the system shows an anomalous 
heavy-fermion behaviour which is accompanied by a large lattice fluctuation and an extreme phonon softening. 
We also calculate an effective potential for the ions and find that a simple harmonic potential for g ^ gc changes 
into a double- well potential for g It gc. The effective pairing interaction between the conduction electrons shows a 
maximum a.t g ~ g^ where the superconducting transition temperature is expected to be maximum. 



The A15 compounds such as VsSi and NbsGe 
have long been attracted much interest as they 
show high Tc and high Hc2 superconductivity as 
well as anomalously large resistivity and Debye- 
Waller factor. Yu and Anderson[l] originally pro- 
posed a local electron-phonon model where the 
strong electron-phonon coupling causes an effec- 
tive double-well potential for the ion which is 
responsible for the anomalous behaviour observed 
in A15 compounds. As a strong coupling fixed 
point, the two- level Kondo systems were investi- 
gated to describe a heavy-fermion like behaviour 
observed in such compounds [2]. More recently, the 
local electron-phonon model have been extensively 
studied by using the NRG approach [3, 4], but pe- 
riodic (lattice) models were not discussed there. 

Recently, another interesting class of materi- 
als has been observed in the filled skutterdites 
such as PrOs4Sbi2[5] and the clathrates such as 
Ce3Pa2oGe6[6], where the rare-earth ion shows a 
rattling motion under a potential with several off- 
center minima. With the new findings, theoretical 
studies on a periodic Anderson model coupled 
with local phonons, i.e., the periodic Anderson- 
Holstein model are highly desirable. The purpose 
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of this paper is to present the results of the dynam- 
ical mean-filed theory (DMFT) for the periodic 
Anderson-Holstein model to elucidate the effect of 
the strong electron-phonon coupling on the heavy- 
fermion behaviour, the effective potential for the 
ions and the superconductivity. 
Our Hamiltonian is given by 

ij(T ia 

ia i 
i i (J 

where c|^, //^ and 6| are creation operators for a 
conduction (c)-electron with spin a at site i, for 
a /-electron and for a phonon, respectively, and 
n/ia = fiafi<y The quantities, e/, V, U and 5, 
are the atomic /-level, the mixing between the c- 
and /-electrons, the on-site Coulomb interaction 
and the electron-phonon coupling strength. The 
density of /-electrons couples with the Einstein 
phonons whose frequency is t^o • 

To solve this model eq.(l), we use the DMFT 
in conjunction with the exact diagonalization 
(ED) method[7]. In the DMFT, the model eq.(l) 
is mapped onto an effective single impurity 
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Anderson- Holstein model[8,9]. Then, the local 
Green's function for a /-electron, 



G/cr(iW„) 



{TrfiAT)flme'"-^dT, (2) 



can be given by the impurity Green's function for 
the impurity Anderson-Holstein model. In the limit 
of infinite dimensions, the self-energy S(ia'„) be- 
comes purely site-diagonal. Then the local Green's 
function for the /-electron satisfies the following 
self-consistency condition: 



de- 



Pie) 



y2 



iw„ - e/ - S(ia;„) 
= [G/<,(ia;„)-i-S(ia;„)]-i, (3) 



where p(e) is the density of states (DOS) for the 
bare conduction band, G / (iw„) is the Green's func- 
tion for the impurity Andcrson-Holstcin model 
with U = g — in an effective medium which will 
be determined self-consistently. In addition, we 
obtain the local Green's function for the phonon. 



.) = -/(r. 



h(r}bU0}y''~'dr. (4) 



To solve the impurity Anderson-Holstein model, 
we use the ED method for a finite-size cluster[9,10] . 
All calculations arc performed at T = 0, and we 
replace the Matsubara frequencies by a fine grid 
of imaginary frequencies LOn = (2n -|- 1)tt//3 with a 
fictitious inverse temperature /? which determines 
the energy resolution. 

The numerical results for 8-site and /3 = 4000 
are almost the same as those for 6-site and (3 = 
200 [11], therefore, we choose the parameters as 8- 
site and /3 = 4000. Also, the numerical results for 
30 phonon states are almost the same as those for 
100 phonon states[8], therefore, we define the cut- 
off of phonon number is 30. In the following numer- 
ical results, we assume a semielliptic DOS for the 
bare conduction band with the bandwidth W = 1, 
p{e) = — e^, and we set wq = 0.05 and V = 
0.2. We concentrate our attention on the particle- 



hole symmetric case with e/ = — where (ric) = 
{nf) = 1. 




Fig. 1. The quasiparticle weight Z, the double-occupancy 
d, the lattice fluctuation (Q^) and the lowest excited energy 
of the phonon spectral function tJo as functions of the 
electron-phonon coupling g for U = (filled circles) and 
U = 0.2 (open circles). {Q^)o and ujo are the corresponding 
values for g = 0. 



1, the quasiparticle weight, Z = 
is plotted as a function of g 



In Fig.l(a 

for {/ = and 0.2. In both cases, Z monoton- 

ically decreases with increasing g and becomes 
extremely small but finite for the strong cou- 



pling g 



where the mass enhancement factor 



m* /m — Z ^ becomes more than one hundred; 
5c ^ 0.085 for f7 = and g^ ^ 0.1 for U = 
0.2 as shown in Fig. 1(a). Thus we can conclude 
that the periodic Anderson-Holstein model eq.(l) 
shows hcavy-fc;rmion behaviour due to the strong 
electron-phonon coupling in the realistic and wide 
parameter range. This is a striking contrast to the 
Holstein-Hubbard model in which the transition 
from the metal with Z ^ {g < Qc) to the bipola- 
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ronic insulator with Z = {g > gc) takes place at 

a critical value of the electron-phonon coupling gc 
and the heavy-fermion behaviour is observed only 
in the narrow range just below gc [8,9] . 

In Fig. 1(b), we plot the double-occupancy, d = 
(n/fn/i), as a function of g ioi U — and 0.2. 
When g increases, d gradually increases for g ^ gc 
while it dose steeply for g ^ gc- The enhancement 
of d is accompanied by the enhancement of the local 
charge fluctuation which causes the heavy fermion 
behaviour observed for g Z gc- 

The lattice fluctuation is defined by {Q^) = 
{Qf — {Qi)"^) with the lattice displacement op- 
erator, Qi = -7= — (bi + b}). Fig. 1(c) shows the 
normalized lattice fluctuation, ((5^)/(Q^)o, where 
(<3^)o = is the lattice fluctuation for g = 
0, i.e., the zero-point oscillation. In the heavy 
fermion regime g ^ gc, we can see the extreme 
enhancement of the lattice fluctuation which is ac- 
companied by the enhancement of the local charge 
fluctuation together with the double occupancy 
shown in Fig. 1(b). 
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Fig. 2. The phonon spectral function — ilmD(ai-|-iO+) for 
several values of 3 at U = 0. 

In Fig. 2, we plot the phonon spectral function, 
-(l/7r)Im£)(u; + iO+), for 5 = 0, 0.06 and 0.08 
at [/ = 0. We find that, with increasing g, the 
multi-phonon state appears and the lowest excited 



energy shifts to low energy. To see the energy shift 

in more detail, we show the lowest excited energy 
uo in the phonon spectral function as a function 
of g in Fig. 1(d). As shown in Fig. 1(d), Qq = ^0 
for g = and ujq decreases with increasing g in 
proportion to Z (see also Fig. 1(a)). A remarkable 
soft phonon mode with « is observed in the 
heavy fermion regime with g ^ gc- 

The large lattice fluctuation and the extreme 
phonon softening observed in Figs. 1(c) and (d) can 
be explained by thinking in terms of an effective 
potential for the ions. To obtain the effective po- 
tential explicitly, we introduce a variational wave 
function for the ions [4], 

|*4Q)|2 = Aexp[-(Bg2 + cq"" + Dq"" + Eq^Jh) 

with q = Q/Qo, where Q is the lattice displace- 
ment and Qo = ,1 is that in the non-interacting 
case {g = 0). In eq.(5), the coefficients A, B, C, 
D and E are the variational parameters which will 
be determined to make (Q2«) = J Q^'^\'^ v{Q)\'^AQ 
{n = 0, 1, 2, 3, 4) as close to (Q^") from the DMFT 
as possible. Then we define the effective potential 
for the ions as 

14ff(0)=log[|vE'„(0)|V|vE'„(Q)|2]. (6) 
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Fig. 3. The effective potential for the ions V^siQ) for 
several values of g at J7 = 0. 
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In Fig.3, Ves{Q) is plotted for several values 
oi g a.t U = 0. In the non-interacting case with 
g = 0, the parameters are A = \/^, B = 0.05 
and C = D = E = which actually yield the 
simple harmonic potential. With increasing g, B 
monotonically decreases while C increases to yield 
an anharmonic potential. In the heavy fermion 
regime with g > g^, B becomes negative while 
keeping C positive to yield a double- well potential 
as shown in Fig.3. This behaviour has already been 
observed in the previous theories [1,3,4] where 
the local electron-phonon model was discussed in 
contrast to the present theory where the periodic 
(lattice) model is discussed. We note that the val- 
ues of D and E are negligibly small and of orders 
of 10-8 - 10-1°. 

Finally, we discuss the superconductivity. The 
renormalized local pairing susceptibility xioc is 
given by 

Xloc = - j {TrC^{T)c^{T)c\{0)c\{Q))dT. (7) 


To calculate Xioc within the ED method, we use 
a spectral representation of r.h.s. in eq.(7) by in- 
serting a complete set of eigenstates[ll]. We also 

calculate the bare local pairing susceptibility Xi'oc 
from the local Green's function Ga{\v) obtained 
from the DMFT, 

X?oc = -rE^THGi(-i^^). (8) 
If 

By using y^^^ and xioc, the local vertex function F 
is given by 

r = [xL]-'-Kioc]-^ (9) 

In the DMFT, the vertex function for the original 
lattice model can be given by the local vertex func- 
tion for the effective impurity model[7]. Therefore, 
r given in eq.(9) corresponds to the effective pair- 
ing interaction for the superconductivity. 

In Fig. 4, we plot the vertex function for /- 
electrons F/ and that for conduction electrons 
Fc together with the quasiparticle weight Z as 
functions of g iov U = 0. When g increases, — F/ 
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Fig. 4. The vertex function for /-electrons T f (filled circles) 
and that for conduction electrons Fc (open circles) together 
with the quasiparticle weight Z (dotted line) as functions 
of g for {/ = 0. 

increases in proportion to the lattice fluctuation 
(Q^) (see also Fig. 1(c)) and dose more rapidly for 
g ^ gc- On the other hand, with increasing g, the 
renormalized band width for /-electrons decreases 
in proportion to Z and becomes almost zero for 
g Z gc- From the both results of — Fj and Z, we 
can suppose the superconducting transition tem- 
perature has a peak around g^ as previously pre- 
dicted by Matsuura and Miyake[2] on the basis of 
the two-level Kondo lattice model. This is directly 
supported by the g-dcpendence of the pairing 
interaction between conduction electrons, — Fc, 
which has a maximum at ~ where the super- 
conducting transition temperature is expected to 
be maximum. 
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